Scar-like structures and their localization in a perfectly square optical billiard 
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We show that scar-hke structures (SLS) in a wide aperture vertical cavity surface emitting laser 
(VCSEL) can be formed even in a perfectly square geometry due to interaction of polarization and 
spatial degrees of freedom of light. We show also that dissipation in the system induces an order 
among the cavity modes, so that SLS become preferred at lasing threshold. More generally, modes 
which are more localized both in coordinate and momentum space have in average lower losses. 
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Quantum billiards traditionally attract a strong at- 
tention P, 0] in the quantum chaos studies. Different 
types of systems, ranging from acoustic and microwave 
resonators to optical cavities and quantum dots belong 
to that class. In most quantum billiards waves of cer- 
tain (not necessary quantum) nature freely move in the 
region of certain shape, surrounded by reflective bound- 
aries. Such systems are described by an energy operator 
H proportional to Laplacian H A = d^x + dyy (in 2D 
case), supplied with corresponding boundary conditions. 

In contrast to billiards, in quantum systems with many 
internal degrees of freedom quantum chaos arises even 
for trivial boundary conditions (if one can speak about 
boundary conditions at all) due to complex structure of 
H. To such class belong nuclei as well as other many- 
body systems [3, Q . We will refer to the later systems 
as to "operator-determined" whereas simple billiards will 
be called "boundary-determined" ones. 

Some systems however belong to an intermediate type, 
for example billiards having a "ball" with nontrivial in- 
ternal structure (i. e. possessing internal degrees of free- 
dom) interacting, in one or another way, with its kinetic 
motion. Up to now only few such systems are known, 
among them are quantum dots in presence of spin-orbital 
coupling [f^, ,5], as well as anisotropic acoustic cavities 
0. In contrast, "photon billiards" (i. e. optical and 
microwave cavities) are traditionally considered as fully 
"boundary-determined" ones flj. 

One important type of an optical billiard is a vertical 
cavity surface emitting laser (VCSEL). Recent advances 
in technology allowed to produce wide aperture, highly 
homogeneous devices of arbitrary shape [3, [1, H, ■ In 
square devices, one of the most prominent features is 
the presence of scar- like structures (SLS) 0, E, 12, 15, 
Ig, |l7|, which are localized along classical trajectories. If 
we consider VCSEL as a "boundary-determined" billiard, 
appearance of such structures must be attributed to some 
other mechanism such as deformations of the boundaries 
Q or mode-locking fld\ . This was done however without 
a rigorous verification. 



provided by a coupling of internal (light polarization) and 
transverse degrees of freedom of photons, which we call 
here polarization-transverse coupling (PTC), appearing 
due to presence of direction-depended anisotropy created 
by the cavity mirrors 

Because fully integrable billiard is separable, the pres- 
ence of SLS points out to a deviation from the complete 
integrability. This deviation can be regulated by mis- 
alignment of the intracavity anisotropy to the bound- 
aries. However, it is always relatively small, making 
the situation similar to quasintegrable scalar billiards 



Moreover, we show that the complicated modes cre- 
ated by PTC are ordered in presence of dissipation, so 
that SLS become preferable (by having less losses) at las- 
ing threshold. More generally, we demonstrate that the 
modes close to threshold are more localized, in average, 
both coordinate and momentum space, comparing to the 
modes with higher losses. This ordering allows to explain 
why SLS so naturally appear in VCSELs. 

Similar relation between localization and dissipation 
was found very recently for a fully chaotic system [isj . 
Localization of long-lived modes appearing in the vicin- 
ity of avoided level crossing (ALC) was also pointed out 
(for the coordinate space only) for dielectric microcavi- 
ties [l^. However, in contrast to [l3|, in our case the 
existence of SLS is neither directly related to ALC, nor 
to dissipation (i. e. to connection to the outer world). 

Using our system as an example we also demonstrate 
that the measures of localization in coordinate and mo- 
mentum space are very sensitive to even a small deviation 
from the complete integrability. 

Despite of sufRciently nonlinear nature of the lasing 
process, many properties of the spatio-temporal distri- 
bution in broad-area VCSELs can be grasped already in 



a linear approximation Although the working area 



In the present article we show that, in fact, the pres- 
ence of SLS does not require any disturbance of the 
square boundaries. Some amount of nonintegrability is 



of VCSEL can be maid very homogeneous in the trans- 
verse direction, the longitudinal structure of the cavity 
is rather complicated. In particular it includes the mul- 
tilayered structure playing the role of the cavity mirrors 
(so called distributed Bragg mirrors, DBRs). However, 
the longitudinal degree of freedom can be excluded from 
the description ^2d\ in an effective way due to single Ion- 
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gitudinal mode operation of the device. As a result of 
such reduction and of subsequent hnearization near the 
lasing threshold the cavity structure is described by 
a single linear operator defining evolution of the complex 
vector field envelope E(rx, t) with time: 



E(ri,t) = 2i/E(ri,i), 



(1) 



where the dot means the partial time derivative, rj^ = 
{x, y} are the transverse coordinates and H is a lin- 
ear operator acting on a transverse field distributions 
E(rj^). H is most easily described for a transversely infi- 
nite VCSEL. In this case, its eigenfunctions are tilted 
waves E ^ g-^r^k^ with certain transverse wavevec- 
tor = {kx,ky}. Therefore, H can be written in 
the transverse Fourier space kj_ = {kx,ky} as a mul- 
tiplication to a 2 X 2 matrix-function /3oo(ki): (3oo = 
ak\ -I- r -t- bs{k±) + iKT{k±). Here a, b and k are 
some constants defined by the parameters of the device, 
k± = |k^p describes a free kinetic motion of the light in- 
side the cavity, F is the intracavity anisotropy, which in 
the Cartesian basis formed by principal anisotropy axis 
can be written as F = diag(7p + i"fa, ^Ip ~ ^Ta), where 
diag(-,-) is a 2 X 2 diagonal matrix with correspond- 
ing elements on the diagonal, 7^ and 7a is the phase 
and amplitude anisotropics. The matrices s(kj^) and 
T(k^) represent the k^-dependent phase and amplitude 



anisotropy, created by DBRs [IJ, l2Qi]. The main axes 
of this anisotropy are perpendicular and parallel to k^ . 
Importantly, T contains the outcoupling losses as well as 
the gain. For the mode at threshold the losses and gain 
exactly compensate each other. 

In the transverse directions the light in VCSEL is 
guided by a thin oxide aperture. Under certain approxi- 
mation (in particular assuming the ideal reflection at the 
side boundaries), the modes of such waveguide are the 
functions of the type E cos{kxx) cos{kyy) which contain 
four spots with equal amplitudes and polarization direc- 
tions in kj^-space. In contrast, eigenmodes of DBRs have 
polarization either perpendicular or parallel to k^ , which 
can not be represented by any combination of waveguide 
modes with a fixed k±. Therefore DBR reflection un- 
avoidly rescatters the eigenmodes of the waveguide into 
the ones with different k±, which creates PTC 12 1. 

Using Poo and the above mentioned properties of the 
waveguide modes, one can directly construct an operator 
which represents H in basis of waveguide modes (see 
[l2l| for details), as well as its reduction (3p with com- 
pletely neglected PTC: 



ijklmn -^ij) , ^ijklmn 



djklmn 



(2) 



(7} 

Here Enm is the jth polarization of the transverse mode 
n, m of the waveguide and Snm is the Kronecker (5-symbol. 

For numerical computation of the eigenvalues and 
eigenfunctions of p]okimn transformed into a square 
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FIG. 1: (Color online) Statistics of the eigenvalues of /3s (red 
solid curve) and /3p (blue dot-dashed curve) for a = 7r/15. 
Dotted green lines shows Poissonian and Wigner distribu- 
tions. In the inset, the statistics P(|Ep) is shown (with 
the same denotations as in the main figure; dotted green line 
shows Porter-Thomas distribution). 



matrix by introducing indices / = i + 2k + 2(iVmax — 
A^min + l)m, J = j + 2l + 2(iV„,ax " iVmin + l)n, whcrc 
Ainax and A'min arc the maximal and minimal waveguide 
mode numbers, defining a cut-off for high- and low- order 
modes. Physically, cut-off for high-order modes is neces- 
sary because they are not guided anymore in transverse 
direction. On the other hand, if wc consider the struc- 
tures formed from the modes which are sufficiently far 
from k±_ — 0, very low order modes can be also neglected. 
For the simulations the values of A^max — 30, A^min = 10 
were taken. In this case the matrix Z?^'"^ has the size 
^ 1000 X 1000. Typical VCSEL parameters were used 
in simulations, in particular, the intracavity anisotropy 
7a = 0.1 ns~^, 7p = 30 ns~^ were taken. The detun- 
ing 6 of the cavity resonance from the gain line center 
(which enters to T) controls the values of k± which have 
maximal gain [2l[. It is chosen large enough so that the 
modes with k± above the high-fc^ cut-off are preferred 
in the infinite device. Due to presence of the cut-off, 
k_L near diagonal \kx\ = \ky\ are selected at threshold 

Sin. 

The statistics P{s) of nearest neighbor separation Si ^ 
Ei+i — Ei of the eigenvalues Ei of the matrices (3s and 
[3p is presented in Fig. [T] in comparison to Poissonian 
(P(s) = e-'') and Wigner (P(s) = f e"'^""/'') distribu- 
tions. The anisotropy axes in Fig. [1] are assumed to be 
rotated to a small angle a — 7r/15 in respect to a;-axis 
(such small displacement may also exist in real devices 
^Mi, although in tendency the anisotropy is aligned to 
the boundaries). In this case, noticeable deviation from 
the Poissonian statistic is present for (3s- On the other 
hand, the eigenvalues of (3p obey Poissonian statistics for 
every a. In contrast, for a = 0, the deviation of P{S) 
from the Poissonian distribution for (3s is not noticeable 
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FIG. 2: (Color online) Example of scared (a)-(e) as well as 
non-localized (f) eigenfunctions of /3s- In (a)-(c) l-Ei^l'^, \Ey'\^ 
and the intensity I of the eigenfunction at threshold (i. e. 
having the lowest losses) are shown; in (d)-(f) the intensity / 
distribution of some subsequent modes is presented. 

anymore ( not shown in Fig. [1]). This shows, that al- 
though PTC itself plays a critical role in the statistics of 
eigenvalues, the alignment of the intracavity anisotropy 
to the boundaries is also important, as it increases the 
degree of mode mixing produced by PTC. 

Another important property of the operator (3^ (both 
for a > and a = 0) is the presence of eigenfunctions 
localized along classical trajectories (i. e. SLS). Exam- 
ples of such structures are shown in Fig. [Ua)-(e). In 
particular, in Fig. [Ua)-(c) the amplitudes of x- and y- 
polarization components as well as the full intensity of 
the mode at threshold (i. e. one with lowest losses) are 
shown. The mode in Fig. ^^c) is localized along tree dif- 
ferent classical trajectories. It is interesting that partic- 
ular polarization components, in contrast to the full in- 
tensity, are not localized along the complete trajectory in 
this case. Similar phenomena were observed also exper- 
imentally [3] (cf- also [ItI where different polarizations 
follow different trajectories). 

In general, SLS are very common for /3s, for both a = 
and a > (see some further examples in Fig. EJd), (e)). 
As a rule, they contain more than one classical trajectory. 
Like [H, [3| , SLS in our system seemingly do not become 
more rear with increasing of fc^. On the other hand, far 
from threshold the eigenfunctions becomes less localized 
(see an example in Fig.[2lf)). 

The existence of SLS clearly points out to deviation 
from the full integrability, because the later assumes an 
existence of coordinate system where the billiard becomes 
fully separable, which excludes the possibility of SLS. For 
a = this deviation is "undetected" by P{s). 

In contrast to (3s, the set of eigenfunctions of (3p do 
not contain SLS. For a = 0, /3p is diagonal. For a ^ it 
consists of 2 X 2 blocks at the main diagonal (describing 
the polarization degrees of freedom). Therefore, all the 
cases we consider can be arranged in order of increasing 
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FIG. 3: (Color online) Localization gk (in momentum space) 

(a) and g (in coordinate space) (b) of the eigenfunctions of fia 
for Of = 7r/15 (red solid line) and for a = (blue dashed line) 
arranged according increasing of the absolute values of the 
imaginary parts of corresponding eigenvalues (i. e. losses). In 

(b) g is also shown for /3p and q = (dotted green line). 



of "nonintegrability" : I3p with a = is the most regular 
(and corresponds to the fully integrable case), whereas 
(3s with a 7^ is the most "nonintegrable" . 

As follows from the previous, the shape of eigenmodes 
can provide a sensitive tool for the description of devi- 
ation of the system from the full integrability. This is 
supported by consideration of the statistics of eigenfunc- 
tion amplitudes P(|Ep) shown in inset to Fig.[T] P(|Ep) 
for (3s for both a 7^ and a = (they are very sim- 
ular to each other and shown by a single red curve in 
the inset) is sufficiently different to the statistics of (3p 
(blue curve). Remarkably, P(|Ep) for (3s is very close to 
Porter-Thomas distribution P(|E|2) = e^l^l'/^/ 
(which is characteristic for chaotic billiards [1|]). 

In this paper however, we are interested in analysis 
not only of different "levels of integrability" , but also of 
the localization prop erties of eigenfunctions. As it was 
shown in 0, [HI, H^, sometimes the good indication of 
SLS is a localization in the "momentum" space : gk — 
AHo/{Ho), where Ho = fci, AHo = ^(-^0) " (^o)^. 
gk becomes exactly zero for the integrable case and is 
shown in Fig.[3l^a) for Ps- One can see that gk is minimal 
at threshold (n ^ 1) in average for a 7^ 0, and grows with 
increasing losses {n — > 00). For "less nonintegrable" case 
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of a = the localization level is nearly constant for low 
n, but starts to grow for n above ~ 200. In general, gk 
can be used to distinguish between SLS and randomly- 
distributed eigenfunctions in chaotic billiards [2^ . but is 
not very well suitable for quasi- integrable billiards, be- 
cause the modes unlocalized in coordinate space can be 
well localized in momentum space [lo| . 

In this letter, we consider also a measure of localiza- 
tion in coordinate space. Namely, we introduce a func- 
tional g[I(x,y)] = J I dxdy / SImax on the spatial inten- 
sity distribution /(x, y) = \Ei{x, y)\^ (here the integration 
is made over the whole billiard area, and the result is nor- 
malized to the area S and to the maximal value of the 
intensity Imax = maxr^ (/(r^)). g varies from one (for 
/ = const) to zero (for a pattern localized near a single 
point, i. e. for the one close to (5-function). Remarkably, 
for an integrable square billiard g = l/A for every eigen- 
function (see Fig. [3l^b), gree dashed curve). Therefore, 
deviation of g from this value may indicate a deviation 
of the system from the complete integrability. In the 
other limit, for fully chaotic case, the typical value of g 
is rather small. One can estimate it using the fact that 
eigenfunctions in fully chaotic billiards can be simulated 
as a sum of the waves with similar fc_L but random phases 
[l|. For such a sum, according to our computations, the 
value of g averaged over large number of eigenfunctions 
depends only on k± and is in the range of 0.06 — 0.07 for 
k±_ used in the present article. 

The values of g for (3^ and (3p are shown in Fig. [31 For 
the case of [3s the deviation of g from the "integrable" 
value 0.25 is quite strong. Intriguingly, in analogy to gu 
the average value of g grows with n. Thus, SLS appear- 
ing at threshold can be described as the patterns having 
lowest g and gu simultaneously, supporting their consid- 
eration as coherent states 11. general, the 
whole pair {g,gk} can be considered as a useful measure, 
allowing to distinguish between SLS, delocalized as well 
as "chaotic" structures. 

It should be noted that the deviation from the full 
integrability in VCSEL disappears in a circular geometry 
, because the terms s and T in /3oo are isotropic in that 
case. In this sense, the situation is similar to the case of 
spin-orbit coupling of electrons in quantum dots 

From the experimental point of view, it is rather prob- 
lematic to distinguish clearly between the effects, appear- 
ing due to small deviations of the boundary conditions 
from a perfect square, and due to deviation of operator 
of the system H from the Laplacian. The result of the 
present article shows that the role of the former is often 
overestimated. 

Whereas localization in the momentum space has a 
concrete physical reason in the case of VCSEL (certain 
kj_ have the highest gain because their frequencies are 
closer to the gain maximum), localization in the coor- 
dinate space is not obvious from the physical point of 
view. Even localization in the momentum space depends 



on the level of integrability of the system (regulated by 
a), although the change of a does not affect the losses 
directly. 

An intriguing similarity to the results of [l^ should 
be pointed out, despite the mechanism of emission in 
VCSEL is completely different from the one in dielectric 
microcavities. In [19|, long lived states created by con- 
nection of the neighboring modes via continuum of ex- 
ternal ones (i. e. mediated by dissipation) in the vicinity 
of ALC were shown to be localized in coordinate space. 
In contrast to [19], appearance of SLS in our system is 
related only to PTC and not to the presence of dissipa- 
tion. Namely, according to our numerical simulations, if 
we remove all the losses in the system (i. e. assume k = 
in /3oo), SLS do not disappear. Moreover, the shape of 
many SLS is not altered significantly. Thus, the presence 
of dissipation does not create SLS but orders them in- 
stead. In addition, SLS are also not directly related to 
ALC in our case, because the frequencies of the modes 
grow in average monotonically with their losses, which is 
not the case in the vicinity of ALC. 

In general, the results of this Letter together with 18|, 
fiol allows us to suspect a general mechanism, which is 
independent from every particular physical realization of 
open billiard, leading to a relation between localization 
and losses, which is still to be clarified. 

The author is grateful for useful discussions with T. 
Ackemann. 
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